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We review the use of Type la supernovae for cosmological distance determinations. Low- 
redshift SNe la (z < 0.1) demonstrate that (a) the Hubble expansion is linear, (b) Ho = 65 ± 2 
(statistical) km s _1 Mpc -1 , (c) the bulk motion of the Local Group is consistent with the COBE 
result, and (d) the properties of dust in other galaxies are similar to those of dust in the Milky 
Way. We find that the light curves of high-redshift (z = 0.3-1) SNe la are stretched in a manner 
consistent with the expansion of space; similarly, their spectra exhibit slower temporal evolution 
(by a factor of 1 + z) than those of nearby SNe la. The luminosity distances of our first set of 
16 high-redshift SNe la are, on average, 10-15% farther than expected in a low mass-density 
(Q.M = 0.2) universe without a cosmological constant. Preliminary analysis of our second set of 
9 SNe la is consistent with this. Our work strongly supports eternally expanding models with 
positive cosmological constant and a current acceleration of the expansion. We address many 
potential sources of systematic error; at present, none of them reconciles the data with Qa = 
and go > 0. The dynamical age of the Universe is estimated to be 14.2 ± 1.7 Gyr, consistent 
with the ages of globular star clusters. 



1. Introduction 

Supernovae (SNe) come in two main varieties (see Filippenko 1997b for a review). 
Those whose optical spectra exhibit hydrogen are classified as Type II, while hydrogen- 
deficient SNe are designated Type I. SNe I are further subdivided according to the ap- 
pearance of the early-time spectrum: SNe la are characterized by strong absorption near 
6150 A (now attributed to Si II), SNe lb lack this feature but instead show prominent 
He I lines, and SNe Ic have neither the Si II nor the He I lines. SNe la are believed to 
result from the thermonuclear disruption of carbon-oxygen white dwarfs, while SNe II 
come from core collapse in massive supergiant stars. The latter mechanism probably 
produces most SNe Ib/Ic as well, but the progenitor stars previously lost their outer 
layers of hydrogen or even helium. 

It has long been recognized that SNe la may be very useful distance indicators for 
a number of reasons (Branch & Tammann 1992; Branch 1998, and references therein). 
(1) They are exceedingly luminous, with peak absolute blue magnitudes averaging —19.2 
if the Hubble constant, Hq, is 65 km s _1 Mpc^" 1 . (2) "Normal" SNe la have small 
dispersion among their peak absolute magnitudes (a <J 0.3 mag). (3) Our understanding 
of the progenitors and explosion mechanism of SNe la is on a reasonably firm physical 
basis. (4) Little cosmic evolution is expected in the peak luminosities of SNe la, and 
it can be modeled. This makes SNe la superior to galaxies as distance indicators. (5) 
One can perform local tests of various possible complications and evolutionary effects by 
comparing nearby SNe la in different environments. 

Research on SNe la in the 1990s has demonstrated their enormous potential as cos- 
mological distance indicators. Although there are subtle effects that must indeed be 
taken into account, it appears that SNe la provide among the most accurate values of 
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Hq, go (the deceleration parameter), Qm (the matter density), and SI a (the cosmological 
constant, Ac 2 /SH§). 

There are now two major teams involved in the systematic investigation of high-rcdshift 
SNe la for cosmological purposes. The "Supernova Cosmology Project" (SCP) is led by 
Saul Pcrlmutter of the Lawrence Berkeley Laboratory, while the "High-Z Supernova 
Search Team" (HZT) is led by Brian Schmidt of the Mt. Stromlo and Siding Springs 
Observatories. One of us (A.V.F.) has worked with both teams, but his primary allegiance 
is now with the HZT. In this lecture we present results from the HZT. 



2. Homogeneity and Heterogeneity 

The traditional way in which SNe la have been used for cosmological distance deter- 
minations has been to assume that they are perfect "standard candles" and to compare 
their observed peak brightness with those of SNe la in galaxies whose distances have 
been independently determined (e.g., Cepheids). The rationale is that SNe la exhibit 
relatively little scatter in their peak blue luminosity (o- B w 0.4-0.5 mag; Branch & Miller 
1993), and even less if "peculiar" or highly reddened objects are eliminated from consid- 
eration by using a color cut. Moreover, the optical spectra of SNe la are usually quite 
homogeneous, if care is taken to compare objects at similar times relative to maximum 
brightness (Riess et al. 1997, and references therein). Branch, Fisher, & Nugent (1993) 
estimate that over 80% of all SNe la discovered thus far are "normal." 

From a Hubble diagram constructed with unreddened, moderately distant SNe la 
(z <; 0.1) for which peculiar motions should be small and relative distances (as given by 
ratios of redshifts) are accurate, Vaughan et al. (1995) find that 

<M B (max)>= (-19.74 ±0.06) + 5 log (iJo/50) mag. (2.1) 

In a series of papers, Sandage et al. (1996) and Saha et al. (1997) combine similar rela- 
tions with Hubble Space Telescope (HST) Ccpheid distances to the host galaxies of seven 
SNe la to derive H = 57 ± 4 km s" 1 Mpc -1 . 

Over the past decade it has become clear, however, that SNe la do not constitute a 
perfectly homogeneous subclass (e.g., Filippenko 1997a, b). In retrospect this should have 
been obvious: the Hubble diagram for SNe la exhibits scatter larger than the photometric 
errors, the dispersion actually rises when reddening corrections are applied (under the 
assumption that all SNe la have uniform, very blue intrinsic colors at maximum; van 
den Bergh & Pazder 1992; Sandage & Tammann 1993), and there are some significant 
outliers whose anomalous magnitudes cannot possibly be explained by extinction alone. 

Spectroscopic and photometric peculiarities have been noted with increasing frequency 
in well-observed SNe la. A striking case is SN 1991T; its pre-maximum spectrum did not 
exhibit Si II or Ca II absorption lines, yet two months past maximum the spectrum was 
nearly indistinguishable from that of a classical SN la (Filippenko et al. 1992b; Phillips 
et al. 1992). The light curves of SN 1991T were slightly broader than the SN la template 
curves, and the object was probably somewhat more luminous than average at maximum. 
The reigning champion of well observed, peculiar SNe la is SN 1991bg (Filippenko et al. 
1992a; Leibundgut et al. 1993; Turatto et al. 1996). At maximum brightness it was 
subluminous by 1.6 mag in V and 2.5 mag in B, its colors were intrinsically red, and 
its spectrum was peculiar (with a deep absorption trough due to Ti II). Moreover, the 
decline from maximum brightness was very steep, the /-band light curve did not exhibit 
a secondary maximum like normal SNe la, and the velocity of the ejecta was unusually 
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low. The photometric heterogeneity among SNe la is well demonstrated by Suntzcff 
(1996) with five objects having excellent BVRI light curves. 



3. Cosmological Uses 

3.1. Luminosity Corrections and Nearby Supernovae 

Although SNe la can no longer be considered perfect "standard candles," they are still 
exceptionally useful for cosmological distance determinations. Excluding those of low 
luminosity (which are hard to find, especially at large distances), most SNe la are nearly 
standard (Branch, Fisher, & Nugent 1993). Also, after many tenuous suggestions (e.g., 
Pskovskii 1977, 1984; Branch 1981), convincing evidence has finally been found for a 
correlation between light-curve shape and luminosity. Phillips (1993) achieved this by 
quantifying the photometric differences among a set of nine well-observed SNe la using a 
parameter, Ami5(B), which measures the total drop (in B magnitudes) from maximum 
to t = 15 days after B maximum. In all cases the host galaxies of his SNe la have 
accurate relative distances from surface brightness fluctuations or from the Tully-Fisher 
relation. In B, the SNe la exhibit a total spread of ~ 2 mag in maximum luminosity, 
and the intrinsically bright SNe la clearly decline more slowly than dim ones. The range 
in absolute magnitude is smaller in V and I, making the correlation with Ami 5 (B) less 
steep than in B, but it is present nonetheless. 

Using SNe la discovered during the Calan/Tololo survey (z < 0.1), Hamuy et al. 
(1995, 1996b) confirm and refine the Phillips (1993) correlation between Amis (2?) and 
M max {B, V): it is not as steep as had been claimed. Apparently the slope is steep only 
at low luminosities; thus, objects such as SN 1991bg skew the slope of the best-fitting 
single straight line. Hamuy et al. reduce the scatter in the Hubble diagram of normal, 
unreddened SNe la to only 0.17 mag in B and 0.14 mag in V; see also Tripp (1997). 

In a similar effort, Riess, Press, & Kirshner (1995a) show that the luminosity of SNe la 
correlates with the detailed shape of the light curve, not just with its initial decline. They 
form a "training set" of light-curve shapes from 9 well-observed SNe la having known 
relative distances, including very peculiar objects (e.g., SN 1991bg). When the light 
curves of an independent sample of 13 SNe la (the Calan/Tololo survey) are analyzed 
with this set of basis vectors, the dispersion in the U-band Hubble diagram drops from 
0.50 to 0.21 mag, and the Hubble constant rises from 53 ± 11 to 67 ± 7 km s -1 Mpc -1 , 
comparable to the conclusions of Hamuy et al. (1995, 1996b). About half of the rise in 
Hq results from a change in the position of the "ridge line" defining the linear Hubble 
relation, and half is from a correction to the luminosity of some of the local calibrators 
which appear to be unusually luminous (e.g., SN 1972E). 

By using light-curve shapes measured through several different filters, Riess, Press, & 
Kirshner (1996a) extend their analysis and objectively eliminate the effects of interstellar 
extinction: a SN la that has an unusually red B — V color at maximum brightness is 
assumed to be intrinsically subluminous if its light curves rise and decline quickly, or of 
normal luminosity but significantly reddened if its light curves rise and decline slowly. 
With a set of 20 SNe la consisting of the Calan/Tololo sample and their own objects, 
Riess, Press, & Kirshner (1996a) show that the dispersion decreases from 0.52 mag to 0.12 
mag after application of this "multi-color light curve shape" (MLCS) method. The results 
from a very recent, expanded set of nearly 50 SNe la indicate that the dispersion decreases 
from 0.44 mag to 0.15 mag (Riess et al. 1999b). The resulting Hubble constant is 65 ± 2 
(statistical) km s -1 Mpc -1 , with an additional systematic and zero-point uncertainty 
of ±5 km s -1 Mpc -1 . Riess, Press, & Kirshner (1996a) also show that the Hubble 
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flow is remarkably linear; indeed, SNe la now constitute the best evidence for linearity. 
Finally, they argue that the dust affecting SNe la is not of circumstellar origin, and show 
quantitatively that the extinction curve in external galaxies typically does not differ from 
that in the Milky Way (cf. Branch & Tammann 1992; but see Tripp 1998). 

Riess, Press, & Kirshncr (1995b) capitalize on another use of SNe la: determination 
of the Milky Way Galaxy's peculiar motion relative to the Hubble flow. They select 
galaxies whose distances were accurately determined from SNe la, and compare their 
observed recession velocities with those expected from the Hubble law alone. The speed 
and direction of the Galaxy's motion are consistent with what is found from COBE 
(Cosmic Background Explorer) studies of the microwave background, but not with the 
results of Lauer & Postman (1994). 

The advantage of systematically correcting the luminosities of SNe la at high redshifts 
rather than trying to isolate "normal" ones seems clear in view of recent evidence that 
the luminosity of SNe la may be a function of stellar population. If the most luminous 
SNe la occur in young stellar populations (Hamuy et al. 1995, 1996a; Branch, Roman- 
ishin, & Baron 1996), then we might expect the mean peak luminosity of high-redshift 
SNe la to differ from that of a local sample. Alternatively, the use of Cepheids (Popula- 
tion I objects) to calibrate local SNe la can lead to a zero point that is too luminous. On 
the other hand, as long as the physics of SNe la is essentially the same in young stellar 
populations locally and at high redshift, we should be able to adopt the luminosity cor- 
rection methods (photometric and spectroscopic) found from detailed studies of nearby 
samples of SNe la. 



4. High-Redshift Supernovae 

4.1. The Search 

These same techniques can be applied to construct a Hubble diagram with high-redshift 
SNe, from which the value of q can be determined. With enough objects spanning a range 
of redshifts, we can determine O jv/ an d independently (e.g., Goobar & Perlmutter 
1995). Contours of peak apparent i?-band magnitude for SNe la at two redshifts have 
different slopes in the Qm-^a plane, and the regions of intersection provide the answers 
we seek. 

Based on the pioneering work of Norgaard-Niclscn et al. (1989), whose goal was to 
find SNe in moderate-redshift clusters of galaxies, Perlmutter et al. (1997) and our HZT 
(Schmidt et al. 1998) devised a strategy that almost guarantees the discovery of many 
faint, distant SNe la on demand, during a predetermined set of nights. This "batch" 
approach to studying distant SNe allows follow-up spectroscopy and photometry to be 
scheduled in advance, resulting in a systematic study not possible with random dis- 
coveries. Most of the searched fields are equatorial, permitting follow-up from both 
hemispheres. 

Our approach is simple in principle; see Schmidt et al. (1998) for details, and for a 
description of our first high-redshift SN la (SN 1995K). Pairs of first-epoch images are 
obtained with the CTIO or CFHT 4-m telescopes and wide-angle imaging cameras during 
the nights just after new moon, followed by second-epoch images 3-4 weeks later. (Pairs 
of images permit removal of cosmic rays, asteroids, and distant Kuiper-belt objects.) 
These are compared immediately using well-tested software, and new SN candidates are 
identified in the second-epoch images. Spectra are obtained as soon as possible after 
discovery to verify that the objects are SNe la and determine their redshifts. Each team 
has already found over 70 SNe in concentrated batches, as reported in numerous IAU 
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Circulars (e.g., Pcrlmuttcr et al. 1995 — 11 SNe with 0.16 £ z £ 0.65; Suntzeff ei a/. 
1996 — 17 SNe with 0.09 ^ z ^ 0.84). 

Intensive photometry of the SNe la commences within a few days after procurement 
of the second-epoch images; it is continued throughout the ensuing and subsequent dark 
runs. In a few cases HST images are obtained. As expected, most of the discoveries 
are on the rise or near maximum brightness. When possible, the SNe are observed in 
filters which closely match the redshifted B and V bands; this way, the if-corrections 
become only a second-order effect (Kim, Goobar, & Perlmutter 1996). Custom-designed 
filters for redshifts centered on 0.35 and 0.45 are used by our HZT (Schmidt et al. 1998), 
when appropriate. We try to obtain excellent multi-color light curves, so that reddening 
and luminosity corrections can be applied (Riess, Press, & Kirshner 1996a; Hamuy et al. 
1996a,b). 

Although SNe in the magnitude range 22-22.5 can sometimes be spectroscopically 
confirmed with 4-m class telescopes, the signal-to-noise ratios are low, even after several 
hours of integration. Certainly Keck is required for the fainter objects (mag 22.5-24.5). 
With Keck, not only can we rapidly confirm a large number of candidate SNe, but we 
can search for peculiarities in the spectra that might indicate evolution of SNe la with 
redshift. Moreover, high-quality spectra allow us to measure the age of a supernova: we 
have developed a method for automatically comparing the spectrum of a SN la with a 
library of spectra corresponding to many different epochs in the development of SNe la 
(Riess et al. 1997). Our technique also has great practical utility at the telescope: we can 
determine the age of a SN "on the fly," within half an hour after obtaining its spectrum. 
This allows us to rapidly decide which SNe are best for subsequent photometric follow-up, 
and we immediately alert our collaborators on other telescopes. 

4.2. Results 

First, we note that the light curves of high-rcdshift SNe la are broader than those of 
nearby SNe la; the initial indications of Leibundgut et al. (1996) and Goldhaber et al. 
(1997) are amply confirmed with our larger samples. Quantitatively, the amount by 
which the light curves are "stretched" is consistent with a factor of 1 + z, as expected 
if redshifts are produced by the expansion of space rather than by "tired light." We 
were also able to demonstrate this spectroscopically at the 2a confidence level for a single 
object: the spectrum of SN 1996bj (z = 0.57) evolved more slowly than those of nearby 
SNe la, by a factor consistent with 1 + z (Riess et al. 1997). More recently, we have used 
observations of SN 1997ex (z = 0.36) at three epochs to conclusively verify the effects of 
time dilation: temporal changes in the spectra are slower than those of nearby SNe la 
by roughly the expected factor of 1.36 (Filippenko et al. 1999). 

Following our Spring 1997 campaign, in which we found a SN with z — 0.97 (SN 
1997ck), and for which we obtained HST follow-up of three SNe, we published our 
first substantial results concerning the density of the Universe (Garnavich et al. 1998a): 
Q.M = 0.35 ± 0.3 under the assumption that fi tota j = 1, or Qm = —0.1 ± 0.5 under the 
assumption that Q\ = 0. Our independent analysis of 10 SNe la using the "snapshot" 
distance method (with which conclusions are drawn from sparsely observed SNe la) gives 
quantitatively similar conclusions (Riess et al. 1998a). 

Our next results, obtained from a total of 16 high-z SNe la, were announced at a confer- 
ence in February 1998 (Filippenko & Riess 1998) and formally published in September 
1998 (Riess et al. 1998b). The Hubble diagram (from a refined version of the MLCS 
method; Riess et al. 1998b) for the 10 best-observed high-z SNe la is given in Figure 
1, while Figure 2 illustrates the derived confidence contours in the I^m-^a plane. We 
confirm our previous suggestion that Qm is low. Even more exciting, however, is our 
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Figure 1. The upper panel shows the Hubble diagram for the low-redshift and high-redshift 
SN la samples with distances measured from the MLCS method; see Riess et al. (1998b). Over- 
plotted are three world models: "low" and "high" Qm with Qa = 0, and the best fit for a flat 
universe (Qm = 0.28, Qa = 0.72). The bottom panel shows the difference between data and 
models from the Qm = 0.20, Qa = prediction. Except for SN 1997ck (open symbol; z = 0.97), 
which lacks spectroscopic confirmation and was excluded from the fit, only the 9 best-observed 
high-redshift SNe la are shown. The average difference between the data and the Qm = 0.20, 
Qa = prediction is 0.25 mag. 

Figure 2. Joint confidence intervals for (Qm,Qa) from SNe la (Riess et al. 1998b). The solid 
contours are results from the MLCS method applied to 10 well-observed SN la light curves, 
together with the snapshot method (Riess et al. 1998a) applied to 6 incomplete SN la light 
curves. The dotted contours are for the same objects excluding SN 1997ck (z = 0.97). Regions 
representing specific cosmological scenarios are illustrated. 

conclusion that Qa is nonzero at the 3a statistical confidence level. With the MLCS 
method applied to the full set of 16 SNe la, our formal results are Qm — 0.24 ± 0.10 if 
Qtotai = 1, or Q.m = —0.35 ± 0.18 (unphysical) if Qa = 0. If wc demand that Qm = 0.2, 
then the best value for Qa is 0.66 ± 0.21. These conclusions do not change significantly 
if only the 9 best-observed SNe la arc used (Fig. 1; Q M = 0.28 ± 0.10 if Q to tai = !)• The 
Atoi5(£?) method yields similar results; if anything, the case for a positive cosmological 
constant strengthens. (For brevity, in this paper we won't quote the Ami5(B) numbers; 
see Riess et al. 1998b.) From an essentially independent set of 42 high-z SNe la (only 2 
objects in common), the SCP obtains almost identical results (Perlmutter et al. 1999). 
This suggests that neither team has made a large, simple blunder! 

Very recently, we have calibrated an additional sample of 9 high-z SNe la, including 
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Figure 3. As in Figure 1, the upper panel shows the Hubble diagram for the low- 2 and high- 2 
SN la samples. Here, we include preliminary analysis of 9 additional SNe la (open squares). 
The bottom panel shows the difference between data and models from the Q.m = 0.20, J7a = 
prediction. 

several observed with HST. Preliminary analysis suggests that the new data are entirely 
consistent with the old results, thereby strengthening their statistical significance. Figure 
3 shows the tentative Hubble diagram; full details will be published elsewhere. 

Though not drawn in Figure 2, the expected confidence contours from measurements 
of the angular scale of the first Doppler peak of the cosmic microwave background radi- 
ation (CMBR) are nearly perpendicular to those provided by SNe la (e.g., Zaldarriaga 
et al. 1997; Eisenstein, Hu, & Tegmark 1998); thus, the two techniques provide com- 
plementary information. The space-based CMBR experiments in the next decade (e.g., 
MAP, Planck) will give very narrow ellipses, but a stunning result is already provided by 
existing measurements (Hancock et al. 1998; Lineweaver & Barbosa 1998): our analysis 
of the data in Riess et al. (1998b) demonstrates that Qm + = 0.94 ± 0.26, when the 
SN and CMBR constraints are combined (Garnavich et al. 1998b; see also Lineweaver 
1998, Efstathiou et al. 1999, and others). As shown in Figure 4, the confidence contours 
are nearly circular, instead of highly eccentric ellipses as in Figure 2. We eagerly look 
forward to future CMBR measurements of even greater precision. 

The dynamical age of the Universe can be calculated from the cosmological parameters. 
In an empty Universe with no cosmological constant, the dynamical age is simply the 
"Hubble time" (i.e., the inverse of the Hubble constant); there is no deceleration. SNe la 
yield H = 65 ± 2 km s _1 Mpc -1 (statistical uncertainty only), and a Hubble time of 
15.1 ± 0.5 Gyr. For a more complex cosmology, integrating the velocity of the expansion 
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Figure 4. The HZT's combined constraints from SNe la (Fig. 1) and the position of the first 
Doppler peak of the CMB angular power spectrum; see Garnavich et al. (1998b). The contours 
mark the 68%, 95.4%, and 99.7% enclosed probability regions. Solid curves correspond to results 
from the MLCS method, while dotted ones are from the Amis(-B) method; all 16 SNe la in 
Riess et al. (1998b) were used. 

from the current epoch (z = 0) to the beginning (z = oo) yields an expression for 
the dynamical age. As shown in detail by Riess et al. (1998b), we obtain a value of 
14. 21^ Gyr using the likely range for ((]m,^a) that we measure. (The precision is 
so high because our experiment is sensitive to roughly the difference between and 
£!a, and the dynamical age also varies approximately this way.) Including the systematic 
uncertainty of the Cepheid distance scale, which may be up to 10%, a reasonable estimate 
of the dynamical age is 14.2 ± 1.7 Gyr. 

This result is consistent with ages determined from various other techniques such as the 
cooling of white dwarfs (Galactic disk > 9.5 Gyr; Oswalt et al. 1996), radioactive dating 
of stars via the thorium and europium abundances (15.2 ± 3.7 Gyr; Cowan et al. 1997), 
and studies of globular clusters (10-15 Gyr, depending on whether Hipparcos parallaxes 
of Cepheids are adopted; Gratton et al. 1997; Chaboyer et al. 1998). Evidently, there is 
no longer a problem that the age of the oldest stars is greater than the dynamical age of 
the Universe. 



5. Discussion 

High-redshift SNe la are observed to be dimmer than expected in an empty Universe 
(i.e., VLm = 0) with no cosmological constant. A cosmological explanation for this obser- 
vation is that a positive vacuum energy density accelerates the expansion. Mass density 
in the Universe exacerbates this problem, requiring even more vacuum energy. For a 
Universe with Q,m = 0.2, the average MLCS distance moduli of the well-observed SNe 
are 0.25 mag larger (i.e., 12.5% greater distances) than the prediction from f^A = 0. 
The average MLCS distance moduli are still 0.18 mag bigger than required for a 68.3% 
(la) consistency for a Universe with Qm — 0.2 and without a cosmological constant. 
The derived value of q is —0.75 ± 0.32, implying that the expansion of the Universe is 
accelerating. It appears that the Universe will expand eternally. 

5.1. Systematic Effects 

A very important point is that the dispersion in the peak luminosities of SNe la (a = 0.15 
mag) is low after application of the MLCS method of Riess et al. (1996a, 1998b). With 16 
SNe la, our effective uncertainty is 0.15/4 w 0.04 mag, less than the expected difference 
of 0.25 mag between universes with Q\ = and 0.76 (and low CIm)', see Figure 1. 
Systematic uncertainties of even 0.05 mag (e.g., in the extinction) are significant, and at 
0.1 mag they dominate any decrease in statistical uncertainty gained with a larger sample 
of SNe la. Thus, our conclusions with only 16 SNe la are already limited by systematic 
uncertainties, not by statistical uncertainties — but of course the 10 new objects further 
strengthen our case. 

Here we explore possible systematic effects that might invalidate our results. Of those 
that can be quantified at the present time, none appears to reconcile the data with 
Oa = 0. Further details can be found in Schmidt et al. (1998) and especially Riess et al. 
(1998b). 
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5.1.1. Evolution 

The local sample of SNe la displays a weak correlation between light-curve shape (or 
luminosity) and host galaxy type, in the sense that the most luminous SNe la with 
the broadest light curves only occur in late-type galaxies. Both early-type and late- 
type galaxies provide hosts for dimmer SNe la with narrower light curves (Hamuy et al. 
1996a). The mean luminosity difference for SNe la in late- type and early- type galaxies 
is ~ 0.3 mag. In addition, the SN la rate per unit luminosity is almost twice as high in 
late-type galaxies as in early-type galaxies at the present epoch (Cappellaro et al. 1997). 
These results may indicate an evolution of SNe la with progenitor age. Possibly relevant 
physical parameters are the mass, metallicity, and C/O ratio of the progenitor (Hoflich, 
Wheeler, & Thielemann 1998). 

We expect that the relation between light-curve shape and luminosity that applies 
to the range of stellar populations and progenitor ages encountered in the late-type and 
early-type hosts in our nearby sample should also be applicable to the range we encounter 
in our distant sample. In fact, the range of age for SN la progenitors in the nearby sample 
is likely to be larger than the change in mean progenitor age over the 4-6 Gyr lookback 
time to the high- z sample. Thus, to first order at least, our local sample should correct 
our distances for progenitor or age effects. 

We can place empirical constraints on the effect that a change in the progenitor age 
would have on our SN la distances by comparing subsamples of low-redshift SNe la 
believed to arise from old and young progenitors. In the nearby sample, the mean 
difference between the distances for the early- type (8 SNe la) and late- type hosts (19 
SNe la), at a given redshift, is 0.04 ± 0.07 mag from the MLCS method. This difference 
is consistent with zero. Even if the SN la progenitors evolved from one population at 
low redshift to the other at high redshift, we still would not explain the surplus in mean 
distance of 0.25 mag over the Oa = prediction. 

Moreover, it is reassuring that initial comparisons of high-redshift SN la spectra ap- 
pear remarkably similar to those observed at low redshift. For example, the spectral 
characteristics of SN 1998ai (z — 0.49) appear to be essentially indistinguishable from 
those of low-redshift SNe la; see Figure 5. In fact, the most obviously discrepant spec- 
trum in this figure is the fourth one from the top, that of SN 1994B (z = 0.09); it is 
intentionally included as a "decoy" that illustrates the degree to which even the spectra 
of nearby SNe la can vary. 

We expect that our local calibration will work well at eliminating any pernicious drift 
in the supernova distances between the local and distant samples. However, we need to 
be vigilant for changes in the properties of SNe la at significant lookback times. Our 
distance measurements could be especially sensitive to changes in the colors of SNe la 
for a given light-curve shape. 

5.1.2. Extinction 

Our SN la distances have the important advantage of including corrections for inter- 
stellar extinction occurring in the host galaxy and the Milky Way. Extinction corrections 
based on the relation between SN la colors and luminosity improve distance precision 
for a sample of nearby SNe la that includes objects with substantial extinction (Riess, 
Press, & Kirshncr 1996a); the scatter in the Hubble diagram is much reduced. Moreover, 
the consistency of the measured Hubble flow from SNe la with late-type and early-type 
hosts (see above) shows that the extinction corrections applied to dusty SNe la at low 
redshift do not alter the expansion rate from its value measured from SNe la in low dust 
environments. 

In practice, our high-rcdshift SNe la appear to suffer negligible extinction; their B — V 
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Figure 5. Spectral comparison (in fx) of SN 1998ai (z = 0.49; Keck spectrum) with low-redshift 
(z < 0.1) SNe la at a similar age (~ 5 days before maximum brightness), from Riess et al. 
(1998b) and Filippenko et al. (1999). The spectra of the low-redshift SNe la were resampled 
and convolved with Gaussian noise to match the quality of the spectrum of SN 1998ai. Overall, 
the agreement in the spectra is excellent, tentatively suggesting that distant SNe la are physically 
similar to nearby SNe la. SN 1994B (z — 0.09) differs the most from the others, and was included 
as a "decoy." 

colors at maximum brightness are normal, suggesting little color excess due to reddening. 
Riess, Press, & Kirshner (1996b) found indications that the Galactic ratios between 
selective absorption and color excess are similar for host galaxies in the nearby (z < 0.1) 
Hubble flow. Yet, what if these ratios changed with lookback time (e.g., Aguirre 1999)? 
Could an evolution in dust grain size descending from ancestral interstellar "pebbles" 
at higher redshifts cause us to underestimate the extinction? Large dust grains would 
not imprint the reddening signature of typical interstellar extinction upon which our 
corrections rely. 

However, viewing our SNe through such grey interstellar grains would also induce a 
dispersion in the derived distances. Using the results of Hatano, Branch, & Deaton 
(1998), Riess et al. (1998b) estimate that the expected dispersion would be 0.40 mag 
if the mean grey extinction were 0.25 mag (the value required to explain the measured 
MLCS distances without a cosmological constant). This is significantly larger than the 
0.21 mag dispersion observed in the high-redshift MLCS distances. Furthermore, most 
of the observed scatter is already consistent with the estimated statistical errors, leaving 
little to be caused by grey extinction. Nevertheless, if we assumed that all of the observed 
scatter were due to grey extinction, the mean shift in the SN la distances would only be 
0.05 mag. With the observations presented here, we cannot rule out this modest amount 
of grey interstellar extinction. 

Grey intergalactic extinction could dim the SNe without either telltale reddening or 
dispersion, if all lines of sight to a given redshift had a similar column density of absorbing 
material. The component of the intergalactic medium with such uniform coverage cor- 



Filippenko & Riess: Cosmological Implications of SNe la 11 

responds to the gas clouds producing Lyman-a forest absorption at low redshifts. These 
clouds have individual H I column densities less than about 10 15 cm~ 2 (Bahcall et al. 

1996) . However, they display low metallicities, typically less than 10% of solar. Grey 
extinction would require larger dust grains which would need a larger mass in heavy 
elements than typical interstellar grain size distributions to achieve a given extinction. 
Furthermore, these clouds reside in hard radiation environments hostile to the survival 
of dust grains. Finally, the existence of grey intergalactic extinction would only augment 
the already surprising excess of galaxies in high-rcdshift galaxy surveys (e.g., Huang et al. 

1997) . 

We conclude that grey extinction does not seem to provide an observationally or phys- 
ically plausible explanation for the observed faintness of high-redshift SNe la. 

5.1.3. Selection Bias 

Sample selection has the potential to distort the comparison of nearby and distant 
SNe. Most of our nearby (z < 0.1) sample of SNe la was gathered from the Calan/Tololo 
survey (Hamuy et al. 1993), which employed the blinking of photographic plates obtained 
at different epochs with Schmidt telescopes, and from less well-defined searches (Riess 
et al. 1999a). Our distant (z ;> 0.16) sample was obtained by subtracting digital CCD 
images at different epochs with the same instrument setup. 

Although selection effects could alter the ratio of intrinsically dim to bright SNe la in 
the nearby and distant samples, our use of the light-curve shape to determine the super- 
nova's luminosity should correct most of this selection bias on our distance estimates. 
Nevertheless, the final dispersion is nonzero, and to investigate its consequences we used 
a Monte Carlo simulation; details are given by Riess et al. (1998b). The results are very 
encouraging, with recovered values of Qm or SIa exceeding the simulated values by only 
0.02-0.03 for these two parameters considered separately. There are two reasons we find 
such a small selection bias in the recovered cosmological parameters. First, the small 
dispersion of our distance indicator (a « 0.15 mag after light-curve shape correction) 
results in only a modest selection bias. Second, both nearby and distant samples include 
an excess of brighter than average SNe, so the difference in their individual selection 
biases remains small. 

Additional work on quantifying the selection criteria of the nearby and distant samples 
is needed. Although the above simulation and others bode well for using SNe la to 
measure cosmological parameters, we must continue to be wary of subtle effects that 
might bias the comparison of SNe la near and far. 

5.1.4. Effect of a Local Void 

Zchavi et al. (1998) find that the SNe la out to 7000 km s _1 may (2-3cr confidence 
level) exhibit an expansion rate which is 6% greater than that measured for the more 
distant objects; see the low-redshift portion of Figure 1. The implication is that the 
volume out to this distance is underdense relative to the global mean density. 

In principle, a local void would increase the expansion rate measured for our low- 
redshift sample relative to the true, global expansion rate. Mistaking this inflated rate 
for the global value would give the false impression of an increase in the low-redshift 
expansion rate relative to the high-rcdshift expansion rate. This outcome could be incor- 
rectly attributed to the influence of a positive cosmological constant. In practice, only 
a small fraction of our nearby sample is within this local void, reducing its effect on the 
determination of the low-redshift expansion rate. 

As a test of the effect of a local void on our constraints for the cosmological parameters, 
we reanalyzed the data discarding the seven SNe la within 7000 km s _1 (d = 108 Mpc 
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for Ho = 65 km s -1 Mpc -1 ). The result was a reduction in the confidence that Qa > 
from 99.7% (3.0<r) to 98.3% (2.4a) for the MLCS method. 

5.1.5. Weak Gravitational Lensing 

The magnification and demagnification of light by large-scale structure can alter the 
observed magnitudes of high-redshift SNe (Kantowski, Vaughan, & Branch 1995). The 
effect of weak gravitational lensing on our analysis has been quantified by Wambsganss 
et al. (1997) and summarized by Schmidt et al. (1998). SN la light will, on average, be 
demagnified by 0.5% at z — 0.5 and 1% at z = 1 in a Universe with a non-negligible 
cosmological constant. Although the sign of the effect is the same as the influence of a 
cosmological constant, the size of the effect is negligible. 

Holz & Wald (1998) have calculated the weak lensing effects on supernova light from 
ordinary matter which is not smoothly distributed in galaxies but rather clumped into 
stars (i.e., dark matter contained in massive compact halo objects). With this scenario, 
microlensing becomes a more important effect, further decreasing the observed supernova 
luminosities at z = 0.5 by 0.02 mag for flM—O-^ (Holz 1998a). Even if most ordinary 
matter were contained in compact objects, this effect would not be large enough to 
reconcile the SNe la distances with the influence of ordinary matter alone. 

5.1.6. Sample Contamination 

Riess et al. (1998b) consider in detail the possibility of sample contamination by SNe 
that arc not SNe la. Of the 16 initial objects, 12 are clearly SNe la and 2 others are 
almost certainly SNe la (though the region near Si II A6150 was poorly observed in the 
latter two). One object (SN 1997ck at z = 0.97) does not have a good spectrum, and 
another (SN 1996E) might be a SN Ic. A reanalysis with only the 14 most probable SNe la 
does not significantly alter our conclusions regarding a positive cosmological constant. 
However, without SN 1997ck we cannot obtain independent values for £Im and f^A- 

Riess et al. (1998b) and Schmidt et al. (1998) discuss several other uncertainties (e.g., 
differences between fitting techniques; if-corrections) and suggest that they, too, are not 
serious in our study. 

5.2. Future Tests 

We are currently trying to increase the sample of high-z SNe la used for measurements 
of cosmological parameters, primarily to quantify systematic effects. Perhaps the most 
probable one is intrinsic evolution of SNe la with redshift: What if SNe la were different 
at high redshift than locally? One way to observationally address this is through careful, 
quantitative spectroscopy of the quality shown in Figure 5 for SN 1998ai. Nearby SNe la 
exhibit a range of spectral properties that correlate with light-curve shape, and we need 
to determine whether high-z SNe la show the same correlations. Comparisons of spectra 
obtained at very early times might be most illuminating: theoretical models predict that 
differences between SNe la should be most pronounced soon after the explosion, when the 
photosphere is in the outermost layers of the ejecta, which potentially undergo different 
amounts of nuclear burning. 

Light curves can also be used to test for evolution of SNe la. For example, one can 
determine whether the rise time (from explosion to maximum brightness) is the same 
for high-z and low-z SNe la; a difference might indicate that the peak luminosities are 
also different (Hoflich, Wheeler, & Thielemann 1998). Moreover, we should see whether 
high-z SNe la show a second maximum in the rest-frame /-band about 25 days after 
the first maximum, as do normal SNe la (e.g., Ford et al. 1993; Suntzeff 1996). Very 
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Figure 6. Fhe HZF SN la data from Figure 1 (open circles) are plotted relative to an empty 
universe (horizontal line). The two faint curves are the best-fitting A model, and the Qm = 1 
(Oa = 0) model. The darker curve shows a systematic bias that increases linearly with z and 
is consistent with our z = 0.5 data. The expected observational uncertainties of hypothetical 
SNe la at redshifts of 0.85 and 1.2 are shown (filled circles). 



subluminous SNe la do not show this second maximum, but rather follow a linear decline 
(Filippenko et al. 1992a). 

One of the most decisive tests to distinguish between A and systematic effects utilizes 
SNe la at z J> 0.85. We observe that z = 0.5 SNe la are fainter than expected even 
in a non-decelerating universe, and interpret this as excessive distances (a consequence 
of positive If so, the deviation of apparent magnitude from the low-f^M, zero- 

A model should actually begin to decrease at z « 0.85; we are looking so far back in 
time that the A effect becomes small compared with Qm (see Fig. 6). If, on the other 
hand, a systematic effect such as non-standard extinction (Aguirre 1998), evolution of the 
white dwarf progenitors (Hoffich, Wheeler, & Thielemann 1998), or gravitational lensing 
(Wambsganss, Cen, & Ostriker 1998) is the culprit, we expect that the deviation of the 
apparent magnitude will continue growing, to first order (Fig. 6). Given the expected 
uncertainties, one or two objects clearly won't be sufficient, but a sample of 10-20 SNe la 
should give a good statistical result. 

Note that a very large sample of SNe la in the range z = 0.3-1.5 could be used 
to see whether A varies with time ("quintessence" models — e.g., Caldwell, Dave, & 
Steinhardt 1998). This is a very ambitious plan, but it should be achievable within 5- 
10 years. Moreover, with a very large sample (200-1000) of high-redshift SNe la, and 
accurate photometry, it should be possible to conduct a "super-MACHO" experiment. 
That is, the light from some SNe la should be strongly amplified by the presence of 
intervening matter, while the majority will be deamplified (e.g., Holz 1998b; Kantowski 
1998). The distribution of amplification factors can be used to determine the type of 
dark matter most prevalent in the Universe (compact objects, or smoothly distributed). 
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6. Conclusions 

When used with care, SNe la are excellent cosmological probes; the current precision 
of individual distance measurements is roughly 5-10%, but a number of subtleties must 
be taken into account to obtain reliable results. SNe la at z <; 0.1 have been used to 
demonstrate that (a) the Hubble flow is definitely linear, (b) the value of the Hubble 
constant is 65 ± 2 (statistical) ±5 (systematic) km s _1 Mpc -1 , (c) the bulk motion of 
the Local Group is consistent with that derived from COBE measurements, and (d) the 
dust in other galaxies is similar to Galactic dust. 

More recently, we have used a sample of 16 high-redshift SNe la (0.16 < z < 0.97) 
to make a number of deductions, as follows. These are supported by our preliminary 
analysis of 9 additional high-z SNe la. 

(1) The luminosity distances exceed the prediction of a low mass-density (Qm ~ 0.2) 
Universe by about 0.25 mag. A cosmological explanation is provided by a positive cos- 
mological constant at the 99.7% (3cr) confidence level, with the prior belief that CIm > 0. 
We also find that the expansion of the Universe is currently accelerating (qo < 0, where 
qo = 51m/2 — ^a), and that the Universe will expand forever. 

(2) The dynamical age of the Universe is 14.2 ±1.7 Gyr, including systematic uncer- 
tainties in the Cepheid distance scale used for the host galaxies of three nearby SNe la. 

(3) These conclusions do not depend on inclusion of SN 1997ck (z = 0.97; uncertain 
classification and extinction) , nor on which of two light-curve fitting methods is used to 
determine the SN la distances. 

(4) The systematic uncertainties presented by grey extinction, sample selection bias, 
evolution, a local void, weak gravitational lensing, and sample contamination currently 
do not provide a convincing substitute for a positive cosmological constant. Further 
studies of these and other possible systematic uncertainties are needed to increase the 
confidence in our results. 

We emphasize that the most recent results of the SCP (Perlmutter et al. 1999) are 
consistent with those of our HZT. This is reassuring — but other, completely independent 
methods are certainly needed to verify these conclusions. The upcoming space CMBR 
experiments hold the most promise in this regard. Although new questions are arising 
(e.g., What is the physical source of the cosmological constant, if nonzero? Are evolving 
cosmic scalar fields a better alternative?), we speculate that this may be the beginning 
of the "end game" in the quest for the values of Am and £1a- 
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